Dual Inhibitors of Main Protease (MPro) and Cathepsin L as Potent Antivirals against SARS-CoV2.

Given the current impact of SARS-CoV2 and COVID-19 on human health and the global economy, the development of direct acting antivirals is of paramount importance. Main protease (MPro), a cysteine protease that cleaves the viral polyprotein, is essential for viral replication. Therefore, MPro is a novel therapeutic target. We identified two novel MPro inhibitors, D-FFRCMKyne and D-FFCitCMKyne, that covalently modify the active site cysteine (C145) and determined cocrystal structures. Medicinal chemistry efforts led to SM141 and SM142, which adopt a unique binding mode within the MPro active site. Notably, these inhibitors do not inhibit the other cysteine protease, papain-like protease (PLPro), involved in the life cycle of SARS-CoV2. SM141 and SM142 block SARS-CoV2 replication in hACE2 expressing A549 cells with IC50 values of 8.2 and 14.7 nM. Detailed studies indicate that these compounds also inhibit cathepsin L (CatL), which cleaves the viral S protein to promote viral entry into host cells. Detailed biochemical, proteomic, and knockdown studies indicate that the antiviral activity of SM141 and SM142 results from the dual inhibition of MPro and CatL. Notably, intranasal and intraperitoneal administration of SM141 and SM142 lead to reduced viral replication, viral loads in the lung, and enhanced survival in SARS-CoV2 infected K18-ACE2 transgenic mice. In total, these data indicate that SM141 and SM142 represent promising scaffolds on which to develop antiviral drugs against SARS-CoV2.

Protein-Ribofuranosyl Interactions Activate Orotidine 5'-Monophosphate Decarboxylase for Catalysis.

The role of a global, substrate-driven, enzyme conformational change in enabling the extraordinarily large rate acceleration for orotidine 5'-monophosphate decarboxylase (OMPDC)-catalyzed decarboxylation of orotidine 5'-monophosphate (OMP) is examined in experiments that focus on the interactions between OMPDC and the ribosyl hydroxyl groups of OMP. The D37 and T100' side chains of OMPDC interact, respectively, with the C-3' and C-2' hydroxyl groups of enzyme-bound OMP. D37G and T100'A substitutions result in 1.4 kcal/mol increases in the activation barrier ΔG⧧ for catalysis of decarboxylation of the phosphodianion-truncated substrate 1-(ß-d-erythrofuranosyl)orotic acid (EO) but result in larger 2.1-2.9 kcal/mol increases in ΔG⧧ for decarboxylation of OMP and for phosphite dianion-activated decarboxylation of EO. This shows that these substitutions reduce transition-state stabilization by the Q215, Y217, and R235 side chains at the dianion binding site. The D37G and T100'A substitutions result in <1.0 kcal/mol increases in ΔG⧧ for activation of OMPDC-catalyzed decarboxylation of the phosphoribofuranosyl-truncated substrate FO by phosphite dianions. Experiments to probe the effect of D37 and T100' substitutions on the kinetic parameters for d-glycerol 3-phosphate and d-erythritol 4-phosphate activators of OMPDC-catalyzed decarboxylation of FO show that ΔG⧧ for sugar phosphate-activated reactions is increased by ca. 2.5 kcal/mol for each -OH interaction eliminated by D37G or T100'A substitutions. We conclude that the interactions between the D37 and T100' side chains and ribosyl or ribosyl-like hydroxyl groups are utilized to activate OMPDC for catalysis of decarboxylation of OMP, EO, and FO.

The role of SERPIN citrullination in thrombosis.

Aberrant protein citrullination is associated with many pathologies; however, the specific effects of this modification remain unknown. We have previously demonstrated that serine protease inhibitors (SERPINs) are highly citrullinated in rheumatoid arthritis (RA) patients. These citrullinated SERPINs include antithrombin, antiplasmin, and t-PAI, which regulate the coagulation and fibrinolysis cascades. Notably, citrullination eliminates their inhibitory activity. Here, we demonstrate that citrullination of antithrombin and t-PAI impairs their binding to their cognate proteases. By contrast, citrullination converts antiplasmin into a substrate. We recapitulate the effects of SERPIN citrullination using in vitro plasma clotting and fibrinolysis assays. Moreover, we show that citrullinated antithrombin and antiplasmin are increased and decreased in a deep vein thrombosis (DVT) model, accounting for how SERPIN citrullination shifts the equilibrium toward thrombus formation. These data provide a direct link between increased citrullination and the risk of thrombosis in autoimmunity and indicate that aberrant SERPIN citrullination promotes pathological thrombus formation.

Crystal Structure of SARS-CoV-2 Main Protease in Complex with the Non-Covalent Inhibitor ML188.

Viral proteases are critical enzymes for the maturation of many human pathogenic viruses and thus are key targets for direct acting antivirals (DAAs). The current viral pandemic caused by SARS-CoV-2 is in dire need of DAAs. The Main protease (Mpro) is the focus of extensive structure-based drug design efforts which are mostly covalent inhibitors targeting the catalytic cysteine. ML188 is a non-covalent inhibitor designed to target SARS-CoV-1 Mpro, and provides an initial scaffold for the creation of effective pan-coronavirus inhibitors. In the current study, we found that ML188 inhibits SARS-CoV-2 Mpro at 2.5 µM, which is more potent than against SAR-CoV-1 Mpro. We determined the crystal structure of ML188 in complex with SARS-CoV-2 Mpro to 2.39 Å resolution. Sharing 96% sequence identity, structural comparison of the two complexes only shows subtle differences. Non-covalent protease inhibitors complement the design of covalent inhibitors against SARS-CoV-2 main protease and are critical initial steps in the design of DAAs to treat CoVID 19.

The Organization of Active Site Side Chains of Glycerol-3-phosphate Dehydrogenase Promotes Efficient Enzyme Catalysis and Rescue of Variant Enzymes.

A comparison of the values of kcat/Km for reduction of dihydroxyacetone phosphate (DHAP) by NADH catalyzed by wild type and K120A/R269A variant glycerol-3-phosphate dehydrogenase from human liver (hlGPDH) shows that the transition state for enzyme-catalyzed hydride transfer is stabilized by 12.0 kcal/mol by interactions with the cationic K120 and R269 side chains. The transition state for the K120A/R269A variant-catalyzed reduction of DHAP is stabilized by 1.0 and 3.8 kcal/mol for reactions in the presence of 1.0 M EtNH3+ and guanidinium cation (Gua+), respectively, and by 7.5 kcal/mol for reactions in the presence of a mixture of each cation at 1.0 M, so that the transition state stabilization by the ternary E·EtNH3+·Gua+ complex is 2.8 kcal/mol greater than the sum of stabilization by the respective binary complexes. This shows that there is cooperativity between the paired activators in transition state stabilization. The effective molarities (EMs) of ∼50 M determined for the K120A and R269A side chains are ≪106 M, the EM for entropically controlled reactions. The unusually efficient rescue of the activity of hlGPDH-catalyzed reactions by the HPi/Gua+ pair and by the Gua+/EtNH3+ activator pair is due to stabilizing interactions between the protein and the activator pieces that organize the K120 and R269 side chains at the active site. This "preorganization" of side chains promotes effective catalysis by hlGPDH and many other enzymes. The role of the highly conserved network of side chains, which include Q295, R269, N270, N205, T264, K204, D260, and K120, in catalysis is discussed.

Human Glycerol 3-Phosphate Dehydrogenase: X-ray Crystal Structures That Guide the Interpretation of Mutagenesis Studies.

Human liver glycerol 3-phosphate dehydrogenase ( hlGPDH) catalyzes the reduction of dihydroxyacetone phosphate (DHAP) to form glycerol 3-phosphate, using the binding energy associated with the nonreacting phosphodianion of the substrate to properly orient the enzyme-substrate complex within the active site. Herein, we report the crystal structures for unliganded, binary E·NAD, and ternary E·NAD·DHAP complexes of wild type hlGPDH, illustrating a new position of DHAP, and probe the kinetics of multiple mutant enzymes with natural and truncated substrates. Mutation of Lys120, which is positioned to donate a proton to the carbonyl of DHAP, results in similar increases in the activation barrier to hlGPDH-catlyzed reduction of DHAP and to phosphite dianion-activated reduction of glycolaldehyde, illustrating that these transition states show similar interactions with the cationic K120 side chain. The K120A mutation results in a 5.3 kcal/mol transition state destabilization, and 3.0 kcal/mol of the lost transition state stabilization is rescued by 1.0 M ethylammonium cation. The 6.5 kcal/mol increase in the activation barrier observed for the D260G mutant hlGPDH-catalyzed reaction represents a 3.5 kcal/mol weakening of transition state stabilization by the K120A side chain and a 3.0 kcal/mol weakening of the interactions with other residues. The interactions, at the enzyme active site, between the K120 side chain and the Q295 and R269 side chains were likewise examined by double-mutant analyses. These results provide strong evidence that the enzyme rate acceleration is due mainly or exclusively to transition state stabilization by electrostatic interactions with polar amino acid side chains.

Enzyme Architecture: Breaking Down the Catalytic Cage that Activates Orotidine 5'-Monophosphate Decarboxylase for Catalysis.

We report the results of a study of the catalytic role of a network of four interacting amino acid side chains at yeast orotidine 5'-monophosphate decarboxylase ( ScOMPDC), by the stepwise replacement of all four side chains. The H-bond, which links the -CH2OH side chain of S154 from the pyrimidine umbrella loop of ScOMPDC to the amide side chain of Q215 in the phosphodianion gripper loop, creates a protein cage for the substrate OMP. The role of this interaction in optimizing transition state stabilization from the dianion gripper side chains Q215, Y217, and R235 was probed by determining the kinetic parameter kcat/ Km for 16 enzyme variants, which include all combinations of single, double, triple, and quadruple S154A, Q215A, Y217F, and R235A mutations. The effects of consecutive Q215A, Y217F, and R235A mutations on Δ G⧧ for wild-type enzyme-catalyzed decarboxylation sum to 11.6 kcal/mol, but to only 7.6 kcal/mol when starting from S154A mutant. This shows that the S154A mutation results in a (11.6-7.6) = 4.0 kcal/mol decrease in transition state stabilization from interactions with Q215, Y217, and R235. Mutant cycles show that ca. 2 kcal/mol of this 4 kcal/mol effect is from the direct interaction between the S154 and Q215 side chains and that ca. 2 kcal/mol is from a tightening in the stabilizing interactions of the Y217 and R235 side chains. The sum of the effects of individual A154S, A215Q, F217Y and A235R substitutions at the quadruple mutant of ScOMPDC to give the corresponding triple mutants, 5.6 kcal/mol, is much smaller than 16.0 kcal/mol, the sum of the effects of the related four substitutions in wild-type ScOMPDC to give the respective single mutants. The small effect of substitutions at the quadruple mutant is consistent with a large entropic cost to holding the flexible loops of ScOMPDC in the active closed conformation.

Primary Deuterium Kinetic Isotope Effects: A Probe for the Origin of the Rate Acceleration for Hydride Transfer Catalyzed by Glycerol-3-Phosphate Dehydrogenase.

Large primary deuterium kinetic isotope effects (1° DKIEs) on enzyme-catalyzed hydride transfer may be observed when the transferred hydride tunnels through the energy barrier. The following 1° DKIEs on kcat/ Km and relative reaction driving force are reported for wild-type and mutant glycerol-3-phosphate dehydrogenase (GPDH)-catalyzed reactions of NADL (L = H, D): wild-type GPDH, ΔΔ G⧧ = 0 kcal/mol, 1° DKIE = 1.5; N270A, 5.6 kcal/mol, 3.1; R269A, 9.1 kcal/mol, 2.8; R269A + 1.0 M guanidine, 2.4 kcal/mol, 2.7; R269A/N270A, 11.5 kcal/mol, 2.4. Similar 1° DKIEs were observed on kcat. The narrow range of 1° DKIEs (2.4-3.1) observed for a 9.1 kcal/mol change in reaction driving force provides strong evidence that these are intrinsic 1° DKIEs on rate-determining hydride transfer. Evidence is presented that the intrinsic DKIE on wild-type GPDH-catalyzed reduction of DHAP lies in this range. A similar range of 1° DKIEs (2.4-2.9) on ( kcat/ KGA, M-1 s-1) was reported for dianion-activated hydride transfer from NADL to glycolaldehyde (GA) [Reyes, A. C.; Amyes, T. L.; Richard, J. P. J. Am. Chem. Soc. 2016, 138, 14526-14529]. These 1° DKIEs are much smaller than those observed for enzyme-catalyzed hydrogen transfer that occurs mainly by quantum mechanical tunneling. These results support the conclusion that the rate acceleration for GPDH-catalyzed reactions is due to the stabilization of the transition state for hydride transfer by interactions with the protein catalyst. The small 1° DKIEs reported for mutant GPDH-catalyzed and for wild-type dianion-activated reactions are inconsistent with a model where the dianion binding energy is utilized in the stabilization of a tunneling ready state.

Orotidine 5'-Monophosphate Decarboxylase: Probing the Limits of the Possible for Enzyme Catalysis.

The mystery associated with catalysis by what were once regarded as protein black boxes, diminished with the X-ray crystallographic determination of the three-dimensional structures of enzyme-substrate complexes. The report that several high-resolution X-ray crystal structures of orotidine 5'-monophosphate decarboxylase (OMPDC) failed to provide a consensus mechanism for enzyme-catalyzed decarboxylation of OMP to form uridine 5'-monophosphate, therefore, provoked a flurry of controversy. This controversy was fueled by the enormous 1023-fold rate acceleration for this enzyme, which had " jolted many biochemists' assumptions about the catalytic potential of enzymes." Our studies on the mechanism of action of OMPDC provide strong evidence that catalysis by this enzyme is not fundamentally different from less proficient catalysts, while highlighting important architectural elements that enable a peak level of performance. Many enzymes undergo substrate-induced protein conformational changes that trap their substrates in solvent occluded protein cages, but the conformational change induced by ligand binding to OMPDC is incredibly complex, as required to enable the development of 22 kcal/mol of stabilizing binding interactions with the phosphodianion and ribosyl substrate fragments of OMP. The binding energy from these fragments is utilized to activate OMPDC for catalysis of decarboxylation at the orotate fragment of OMP, through the creation of a tight, catalytically active, protein cage from the floppy, open, unliganded form of OMPDC. Such utilization of binding energy for ligand-driven conformational changes provides a general mechanism to obtain specificity in transition state binding. The rate enhancement that results from the binding of carbon acid substrates to enzymes is partly due to a reduction in the carbon acid p Ka that is associated with ligand binding. The binding of UMP to OMPDC results in an unusually large >12 unit decrease in the p Ka = 29 for abstraction of the C-6 substrate hydrogen, due to stabilization of an enzyme-bound vinyl carbanion, which is also an intermediate of OMPDC-catalyzed decarboxylation. The protein-ligand interactions operate to stabilize the vinyl carbanion at the enzyme active site compared to aqueous solution, rather than to stabilize the transition state for the concerted electrophilic displacement of CO2 by H+ that avoids formation of this reaction intermediate. There is evidence that OMPDC induces strain into the bound substrate. The interaction between the amide side chain of Gln-215 from the phosphodianion gripper loop and the hydroxymethylene side chain of Ser-154 from the pyrimidine umbrella of ScOMPDC position the amide side chain to interact with the phosphodianion of OMP. There are no direct stabilizing interactions between dianion gripper protein side chains Gln-215, Tyr-217, and Arg-235 and the pyrimidine ring at the decarboxylation transition state. Rather these side chains function solely to hold OMPDC in the catalytically active closed conformation. The hydrophobic side chains that line the active site of OMPDC in the region of the departing CO2 product may function to stabilize the decarboxylation transition state by providing hydrophobic solvation of this product.

Enzyme Architecture: The Role of a Flexible Loop in Activation of Glycerol-3-phosphate Dehydrogenase for Catalysis of Hydride Transfer.

The side chain of Q295 of glycerol-3-phosphate dehydrogenase from human liver ( hlGPDH) lies in a flexible loop, that folds over the phosphodianion of substrate dihydroxyacetone phosphate (DHAP). Q295 interacts with the side-chain cation from R269, which is ion-paired to the substrate phosphodianion. Kinetic parameters kcat/ Km (M-1 s-1) and kcat/ KGA KHPi (M-2 s-1) were determined, respectively, for catalysis of the reduction of DHAP and for dianion activation of catalysis of reduction of glycolaldehyde (GA) catalyzed by wild-type, Q295G, Q295S, Q295A, and Q295N mutants of hlGPDH. These mutations result in up to a 150-fold decrease in ( kcat/ Km)DHAP and up to a 2.7 kcal/mol decrease in the intrinsic phosphodianion binding energy. The data define a linear correlation with slope 1.1, between the intrinsic phosphodianion binding energy and the intrinsic phosphite dianion binding energy for activation of hlGPDH-catalyzed reduction of GA, that demonstrates a role for Q295 in optimizing this dianion binding energy. The R269A mutation of wild-type GPDH results in a 9.1 kcal/mol destabilization of the transition state for reduction of DHAP, but the same R269A mutation of N270A and Q295A mutants result in smaller 5.9 and 4.9 kcal/mol transition-state destabilization. Similarly, the N270A or Q295A mutations of R269A GPDH each result in large falloffs in the efficiency of rescue of the R269A mutant by guanidine cation. We conclude that N270, which interacts for the substrate phosphodianion and Q295, which interacts with the guanidine side chain of R269, function to optimize the apparent transition-state stabilization provided by the cationic side chain of R269.

Enzyme Architecture: Erection of Active Orotidine 5'-Monophosphate Decarboxylase by Substrate-Induced Conformational Changes.

Orotidine 5'-monophosphate decarboxylase (OMPDC) catalyzes the decarboxylation of 5-fluoroorotate (FO) with kcat/Km = 1.4 × 10-7 M-1 s-1. Combining this and related kinetic parameters shows that the 31 kcal/mol stabilization of the transition state for decarboxylation of OMP provided by OMPDC represents the sum of 11.8 and 10.6 kcal/mol stabilization by the substrate phosphodianion and the ribosyl ring, respectively, and an 8.6 kcal/mol stabilization from the orotate ring. The transition state for OMPDC-catalyzed decarboxylation of FO is stabilized by 5.2, 7.2, and 9.0 kcal/mol, respectively, by 1.0 M phosphite dianion, d-glycerol 3-phosphate and d-erythritol 4-phosphate. The stabilization is due to the utilization of binding interactions of the substrate fragments to drive an enzyme conformational change, which locks the orotate ring of the whole substrate, or the substrate pieces in a caged complex. We propose that enzyme-activation is a possible, and perhaps probable, consequence of any substrate-induced enzyme conformational change.

A reevaluation of the origin of the rate acceleration for enzyme-catalyzed hydride transfer.

There is no consensus of opinion on the origin of the large rate accelerations observed for enzyme-catalyzed hydride transfer. The interpretation of recent results from studies on hydride transfer reactions catalyzed by alcohol dehydrogenase (ADH) focus on the proposal that the effective barrier height is reduced by quantum-mechanical tunneling through the energy barrier. This interpretation contrasts sharply with the notion that enzymatic rate accelerations are obtained through direct stabilization of the transition state for the nonenzymatic reaction in water. The binding energy of the dianion of substrate DHAP provides 11 kcal mol-1 stabilization of the transition state for the hydride transfer reaction catalyzed by glycerol-3-phosphate dehydrogenase (GPDH). We summarize evidence that the binding interactions between (GPDH) and dianion activators are utilized directly for stabilization of the transition state for enzyme-catalyzed hydride transfer. The possibility is considered, and then discounted, that these dianion binding interactions are utilized for the stabilization of a tunnel ready state (TRS) that enables efficient tunneling of the transferred hydride through the energy barrier, and underneath the energy maximum for the transition state. It is noted that the evidence to support the existence of a tunnel-ready state for the hydride transfer reactions catalyzed by ADH is ambiguous. We propose that the rate acceleration for ADH is due to the utilization of the binding energy of the cofactor NAD+/NADH in the stabilization of the transition state for enzyme-catalyzed hydride transfer.

Structure-Reactivity Effects on Intrinsic Primary Kinetic Isotope Effects for Hydride Transfer Catalyzed by Glycerol-3-phosphate Dehydrogenase.

Primary deuterium kinetic isotope effects (1°DKIE) on (kcat/KGA, M-1 s-1) for dianion (X2-) activated hydride transfer from NADL to glycolaldehyde (GA) catalyzed by glycerol-3-phosphate dehydrogenase were determined over a 2100-fold range of enzyme reactivity: (X2-, 1°DKIE); FPO32-, 2.8 ± 0.1; HPO32-, 2.5 ± 0.1; SO42-, 2.8 ± 0.2; HOPO32-, 2.5 ± 0.1; S2O32-, 2.9 ± 0.1; unactivated; 2.4 ± 0.2. Similar 1°DKIEs were determined for kcat. The observed 1°DKIEs are essentially independent of changes in enzyme reactivity with changing dianion activator. The results are consistent with (i) fast and reversible ligand binding; (ii) the conclusion that the observed 1°DKIEs are equal to the intrinsic 1°DKIE on hydride transfer from NADL to GA; (iii) similar intrinsic 1°DKIEs on GPDH-catalyzed reduction of the substrate pieces and the whole physiological substrate dihydroxyacetone phosphate. The ground-state binding interactions for different X2- are similar, but there are large differences in the transition state interactions for different X2-. The changes in transition state binding interactions are expressed as changes in kcat and are proposed to represent changes in stabilization of the active closed form of GPDH. The 1°DKIEs are much smaller than observed for enzyme-catalyzed hydrogen transfer that occurs mainly by quantum-mechanical tunneling.

Enzyme Architecture: Self-Assembly of Enzyme and Substrate Pieces of Glycerol-3-Phosphate Dehydrogenase into a Robust Catalyst of Hydride Transfer.

The stabilization of the transition state for hlGPDH-catalyzed reduction of DHAP due to the action of the phosphodianion of DHAP and the cationic side chain of R269 is between 12.4 and 17 kcal/mol. The R269A mutation of glycerol-3-phosphate dehydrogenase (hlGPDH) results in a 9.1 kcal/mol destabilization of the transition state for enzyme-catalyzed reduction of dihydroxyacetone phosphate (DHAP) by NADH, and there is a 6.7 kcal/mol stabilization of this transition state by 1.0 M guanidine cation (Gua+) [J. Am. Chem. Soc. 2015, 137, 5312-5315]. The R269A mutant shows no detectable activity toward reduction of glycolaldehyde (GA), or activation of this reaction by 30 mM HPO32-. We report the unprecedented self-assembly of R269A hlGPDH, dianions (X2- = FPO32-, HPO32-, or SO42-), Gua+ and GA into a functioning catalyst of the reduction of GA, and fourth-order reaction rate constants kcat/KGAKXKGua. The linear logarithmic correlation (slope = 1.0) between values of kcat/KGAKX for dianion activation of wildtype hlGPDH-catalyzed reduction of GA and kcat/KGAKXKGua shows that the electrostatic interaction between exogenous dianions and the side chain of R269 is not significantly perturbed by cutting hlGPDH into R269A and Gua+ pieces. The advantage for connection of hlGPDH (R269A mutant + Gua+) and substrate pieces (GA + HPi) pieces, (ΔGS‡)HPi+E+Gua = 5.6 kcal/mol, is nearly equal to the sum of the advantage to connection of the substrate pieces, (ΔGS‡)GA+HPi = 3.3 kcal/mol, for wildtype hlGPDH-catalyzed reaction of GA + HPi, and for connection of the enzyme pieces, (ΔGS‡)E+Gua = 2.4 kcal/mol, for Gua+ activation of the R269A hlGPDH-catalyzed reaction of DHAP.

Enzyme Architecture: A Startling Role for Asn270 in Glycerol 3-Phosphate Dehydrogenase-Catalyzed Hydride Transfer.

The side chains of R269 and N270 interact with the phosphodianion of dihydroxyacetone phosphate (DHAP) bound to glycerol 3-phosphate dehydrogenase (GPDH). The R269A, N270A, and R269A/N270A mutations of GPDH result in 9.1, 5.6, and 11.5 kcal/mol destabilization, respectively, of the transition state for GPDH-catalyzed reduction of DHAP by the reduced form of nicotinamide adenine dinucleotide. The N270A mutation results in a 7.7 kcal/mol decrease in the intrinsic phosphodianion binding energy, which is larger than the 5.6 kcal/mol effect of the mutation on the stability of the transition state for reduction of DHAP; a 2.2 kcal/mol stabilization of the transition state for unactivated hydride transfer to the truncated substrate glycolaldehyde (GA); and a change in the effect of phosphite dianion on GPDH-catalyzed reduction of GA, from strongly activating to inhibiting. The N270A mutation breaks the network of hydrogen bonding side chains, Asn270, Thr264, Asn205, Lys204, Asp260, and Lys120, which connect the dianion activation and catalytic sites of GPDH. We propose that this disruption dramatically alters the performance of GPDH at these sites.

Enzyme architecture: optimization of transition state stabilization from a cation-phosphodianion pair.

The side chain cation of R269 lies at the surface of l-glycerol 3-phosphate dehydrogenase (GPDH) and forms an ion pair to the phosphodianion of substrate dihydroxyacetone phosphate (DHAP), which is buried at the nonpolar protein interior. The R269A mutation of GPDH results in a 110-fold increase in K(m) (2.8 kcal/mol effect) and a 41,000-fold decrease in k(cat) (6.3 kcal/mol effect), which corresponds to a 9.1 kcal/mol destabilization of the transition state for GPDH-catalyzed reduction of DHAP by NADH. There is a 6.7 kcal/mol stabilization of the transition state for the R269A mutant GPDH-catalyzed reaction by 1.0 M guanidinium ion, and the transition state for the reaction of the substrate pieces is stabilized by an additional 2.4 kcal/mol by their covalent attachment at wildtype GPDH. These results provide strong support for the proposal that GPDH invests the 11 kcal/mol intrinsic phosphodianion binding energy of DHAP in trapping the substrate at a nonpolar active site, where strong electrostatic interactions are favored, and obtains a 9 kcal/mol return from stabilizing interactions between the side chain cation and transition state trianion. We propose a wide propagation for the catalytic motif examined in this work, which enables strong transition state stabilization from enzyme-phosphodianion pairs.

The activating oxydianion binding domain for enzyme-catalyzed proton transfer, hydride transfer, and decarboxylation: specificity and enzyme architecture.

The kinetic parameters for activation of yeast triosephosphate isomerase (ScTIM), yeast orotidine monophosphate decarboxylase (ScOMPDC), and human liver glycerol 3-phosphate dehydrogenase (hlGPDH) for catalysis of reactions of their respective phosphodianion truncated substrates are reported for the following oxydianions: HPO3(2-), FPO3(2-), S2O3(2-), SO4(2-) and HOPO3(2-). Oxydianions bind weakly to these unliganded enzymes and tightly to the transition state complex (E·S(‡)), with intrinsic oxydianion Gibbs binding free energies that range from -8.4 kcal/mol for activation of hlGPDH-catalyzed reduction of glycolaldehyde by FPO3(2-) to -3.0 kcal/mol for activation of ScOMPDC-catalyzed decarboxylation of 1-ß-d-erythrofuranosyl)orotic acid by HOPO3(2-). Small differences in the specificity of the different oxydianion binding domains are observed. We propose that the large -8.4 kcal/mol and small -3.8 kcal/mol intrinsic oxydianion binding energy for activation of hlGPDH by FPO3(2-) and S2O3(2-), respectively, compared with activation of ScTIM and ScOMPDC reflect stabilizing and destabilizing interactions between the oxydianion -F and -S with the cationic side chain of R269 for hlGPDH. These results are consistent with a cryptic function for the similarly structured oxydianion binding domains of ScTIM, ScOMPDC and hlGPDH. Each enzyme utilizes the interactions with tetrahedral inorganic oxydianions to drive a conformational change that locks the substrate in a caged Michaelis complex that provides optimal stabilization of the different enzymatic transition states. The observation of dianion activation by stabilization of active caged Michaelis complexes may be generalized to the many other enzymes that utilize substrate binding energy to drive changes in enzyme conformation, which induce tight substrate fits.